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Abstract 

We give a solution to the inverse moment problem for a certain class of Hessenberg and 
symmetric matrices related to integrable lattices of Toda type. 

1 Introduction 

This paper continues our previous work , where we dealt with the family of integrable Hamil- 
tonian systems in M^" parametrized by index sets I = {ii, . . . ,ik : 1 < ii < . . . < ik = n} and 
generated by Hamiltonians 

n, n k—1 

Hj[Q,p) = -Y^pf + P,e«»+i-«' + ^6^5'^+^-^^ . (1.1) 

i=l l<.i<n;i^ii,... j=l 

This family contains (after an appropriate coordinate changes) such important integrable systems 
as the standard and relativistic Toda lattices, Volterra lattice and lattices of the Ablowitz- Ladik 
hierarchy. 

In we argued that the full Kostant-Toda flows on Hessenberg matrices provide a convenient 
framework to study systems generated by (1.1). Namely, each of these systems possesses a Lax 



representation with a Lax operator given by an n x n upper Hessenberg matrix X/ = Xi{Q, P) that 
belongs to a certain (2n — 2)-dimensional co-adjoint orbit of the upper triangular group. This orbit 
is determined by / and the value of Tr{Xi) = Y^^=i ^i- Furthermore, Tr{Xj) = Hj and Hamilton 
equations of motion are equivalent to the Toda flow Xj = [Xj, {Xi)<o], where (^)<o denotes the 
lower triangular part of a matrix A. The orbit contains a dense open set of elements that admit 
a factorization into a product of elementary bi-diagonal matrices. Written in terms of parameters 
of this factorization, equations of motion become a particular case of the constrained KP lattice 



studied in |13|. 



Each of the systems described above can be linearized via the Moser map X — > m(A, X) = 
((Al - X)-iei, ei) = E°^o ^ where ei = (1, 0, . . . , 0) and sjiX) = (X^e^ei) (see 0, |, |, |). 
Moreover, as we have shown in |^ the Moser map is very useful in establishing a multi-Hamiltonian 
structure of these systems. However, explicit formulas for the inverse of the Moser map seem 
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to be known only in two cases. If Xj is tri-diagonal (/ = {2, . . . , n}) they give a solution in 
terms of Hankel determinants to the classical finite-dimensional moment problem (see, e.g. IH). 
In the "opposite" case that corresponds to the relativistic Toda lattice ( I = {n}) the solution 
to the inverse problem is given in terms of Toeplitz determinants constructed from the moments 
Sj{Xi),j = -n, ... ,n (see 

The main purpose of this paper is to give a solution to the inverse problem for any /. Explicit 
formulas (given, once again, in terms of Toeplitz determinants) are obtained in sect. 3. 

As in the Hessenberg case, the inverse problem for elementary co-adjoint orbits in the symmetric 
case was previously studied for / = {2, . . . , n} and / = {n}. The latter case corresponds to peakons 
solutions of the shallow water equation and was recently comprehensively studied in We treat 
the symmetric case in sect. 4 for an arbitrary /. 

We would like to express our gratitude to Yu. Suris for valuable comments and suggestions. 

2 Elementary Orbits in the Hessenberg Case 

Everywhere below we denote by ejk an elementary matrix (61^6^)^ ^^-^ and by Cj a column vector 

iSa)a=i of the standard basis in M". 

Denote by J an n x n matrix with Is on the first sub-diagonal and Os everywhere else. Let 
b+, n+, b_, n_ be, resp., algebras of upper triangular, strictly upper triangular, lower triangular and 
strictly lower triangular matrices. Denote by TC the set J -H b+ of upper Hessenberg matrices. 

For any matrix A we write its decomposition into a sum of lower triangular and strictly upper 
triangular matrices as 

A = A_+Ao + A+ 

and define A>o = Aq + A+,A<o = Aq + A., Agym = A+ + Aq + A^. 

A linear Poisson structure on 7i is obtained as a pull-back of the Kirillov-Kostant structure 
on b'^, the dual of b_, if one identifies b'^ and 7i via the trace form. A Poisson bracket of two 
functions /i , /2 on then reads 

{/i,/2}(X) = (X,[(V/i(X))<o,(V/2(X))<o]) , (2.1) 

where we denote by {X, Y) the trace form Tiace{XY) and gradients are computed w.r.t. this form. 
Symplectic leaves of the bracket ( |2.1| ) are orbits of the coadjoint action of the group B_ of lower 
triangular invertible matrices: 

Oxo = {J + (Ad„Xo)>o : n e B_} . (2.2) 

Following [P, consider a family of orbits whose members are parameterized by increasing se- 
quences of natural numbers I = {ii, . . . ,ik : 1 < ii < . . . < ik = n}. To each sequence / there 
corresponds a 1-parameter family of 2(n — l)-dimensional coadjoint orbits 

Mj = U^^kOxj+u1 C n , (2.3) 

where 

fc-i 

Xi = eiii + a-i^^, + J . (2.4) 



2 



We call orbits Oxj+ui elementary. 
The set M| of elements of the form 

X = {J + D){1 - Ck)-\l - Ck-i)-^ •••(!- Ci)-\ (2.5) 

where D = diag((ii, ... ,dn) 

ij-i 

Cj = ^ Caea,a+1 , (2.6) 
a=ij_i 

is dense in Mj. 

The following formulae express entries xim {I < m) of X in terms of Cj, di (see f^) : 



a^Zm — dlUlm + Ul-l,m 



{di + Q_i)q • • • Cm-i, <l <m<ij 

dij_iCij_^ • • • Cm— 1) ^j — 1 — I <^ Tn ^ ij 

di+ Q-i, l=m 

0, otherwise 



(Here cq = 0.) 

Define a sequence ei, . . . e„ by setting 



if i = ij for some < j < k 

1 otherwise 



(2.7) 



(2. 



Then X can also be written as[| 

X = (J + L>)(l + ?7i)(l-C/2)-\ (2.9) 

where 

71—1 n—1 
Ui = ^(1 - ea)Caea,a+l > ^2 = ^ eaCaCa^a+i (2.10) 
a=l a=l 

In what follows, we will also use a sequence of integers (i^i)r=i defined by 

i-l 

u,=i{l-e,)-^ef3 (2.11) 
/3=1 

It is easy to check that 

j if i = ij for some < j < k 

- 6/3 otherwise 

It follows from ( ^l2|) that 

(i) sequences (ej) and (ui) uniquely determine each other. 

(ii) If we define, for every i G {I,-- - ,n}, a set Ni = {v^ : q = 1,... ,i} and a number 
ki = max{j : ij < i} = i — X]/3=i ^/S; then 

13=1 13=1 



(2.12) 



^This was suggested to us by Yu. Suris |14[ 
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3 Solution of the Inverse Problem 

Assume that all di ^ and then define the moment sequence S = {si,i G oi X: 

Si = Si{X) = eJX'ei. (3.1) 

Our goal is to express coefficients Ci,di in terms of S. In fact, only a segment Sk+i-n, ■ ■ ■ , Sn+k 
will be needed. 

Let p{X) = ipo{X) := l,pi(A), . . . ,pn-i{X)) be a solution of the truncated eigenvalue problem 

i 

{p{X)X)i = ^ XaiPa~iW + Pi{X) = Xpi-i{X), i = 1, . . . , n - 1 (3.2) 

a=l 

or, equivalently, 

{p{X){J + D){1 + C/i)) . = A ((1 - U2)p), (A), i = 1, . . . , n - 1 . (3.3) 

Clearly, such solution exists for every A and is uniquely defined. Moreover, each pi{X) is a monic 
polynomial of degree i. We can re- write ( |3.3D as a 3-term recursion for polynomials Pi{X) : 

Pi+i{X) + bi+iPi{X) + (1 - ei)aiPi^i{X) = X{pi{X) - eiCiPi-i{X), ) i = 0, . . . ,n - 1 , (3.4) 

where 

bi = di + {I - ei_i)ci_i, ai = diCi . (3.5) 

Lemma 3.1 For any upper Hessenberg matrix X, polynomials Pi{X) defined by jg.^j satisfy 

Pj(X)ei = Cj+i, i = 0, . . . ,n - 1. (3.6) 



Proof. Since 

i 

XCi — ^ ^ Xiy,i(ia ~\~ 
a=l 

and, by (P), 

i 

Xpi^,{X) = Y^ {X)ei + pi{X)ei , 

a=l 

one concludes that sequences of vectors (ei)"^-^ and {pi~i{X)ei)^^-^ are defined by the same recur- 
rence relations. □ 

Lemma 3.2 For i = 1,... ,n — 1, a subspace generated by vectors {e'^Y^^i coincides with a 
subspace generated by vectors (ef X'^")^^;^, where are defined in ( \2.11\ ). In other words, for 
some constants 7i ^ 

^ j efX^ if i = ij 

^^^^ = I efX-Sri.^. otherw^se ^^''^ ^^''^ 

Moreover, for i = 2, . . . ,n, 

(di ■ ■ ■ di^^ 
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Proof. For any X given by ( |2.5| ), ( |2.g| ) , consider an upper triangular matrix 

V = Dil- Ck)-\1 - Ck-iY^ •••(!- Ci)-i = D{1 + Ui){l - U2)-' . 

(If D is invertible, then V is an upper triangular factor in the Gauss factorization of X." 
By (1^), efCj = for / < ij^i and l>ij. Thus, for j = 0, . . . , k, 

e[,_,V = d.,_,e5_^(l - •••(!- Ci)-^ = d,^_,ef^_^{l - C,)-i(mod 

= di^_^Ci^_, ■ ■ ■ Cj^-ieJ (mod . 

Similar argument shows that efV& -^ij-i ^'^'^ ^ ^ ^j-i- This implies 



e(V^ = e{V^ (mod -Cj.-i) = ( n;0=o ^i/s ' ' -Ci^+i-i) % (mod Z 



cidi nK ^/s^"") 4 (3.9) 

On the other hand, for / G {ij-i, ■ ■ ■ i^j ~ 1}; define m > to be a number such that / + /? = 

for P = 1, . . . ,m and / + m + 1 < ij+m- In other words, / + m + 1 is the smallest index greater than 

I that does not belong to the index set /. Then 

efV-' = efil - C,) ■■■{!- Ck)D-^ = {{-ir+'ci ■ ■ ■ ci+^d-^„^^^) ef+„+i(mod . (3.10) 

Let us denote by J = {/i < . . . < Zn-fc-i} the set {!,... ,n}\I. Then ( 3.10| ) implies 

efX-° = elV-" (mod 
'(_iy.~ici...Q_id,;i---dr,')€ (mod£i_i) . (3.11) 



Note now that = i if and only if a = e^. Furthermore, d^^"*^ • • • di^ = \Yp=2 dp'^'^ ■ Thus, one 

can re-write ( |3.11| ) as 



^ V=2C,dl- '^ 

di---di 



efy-^ri^- = (-1)^-^ n ^Pd7+Tel (mod £._i) = (-l)^-^cidi%^^^ (n,od 

13=1 



(3.12) 



Combining (U), (|1|) with (|2T^ ), one concludes that 

elX""^ = 7ief (mod £^^1) , (3.13) 



where ji is defined by (3.8). This implies the statement of the lemma. □ 
Corollary 3.3 

i i 

eJxy,{X)ei = 0, 1 - ^ < a < i - ^ e;3 (3.14) 
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Proof. By ( |3?7D, e ^X'^'e^+i = for Z = 1, . . . , i. But, by Lemma 2.1, a+i = pi{X)ei. Then (|3l4D 
follows from ( |2.13| ). □ 
Define Toeplitz matrices 



and Toeplitz determinants 



A 



.(0 



det T, 



(s/+a-/3)a,/3=l 



(3.15) 



(3.16) 



Let us also define polynomials 



Sl Sl+l ■ ■ ■ Sl+i 

Sl-i+l Sl-i+2 ■ ■ ■ 

1 A •••A' 



(3.17) 



V^''^ is a determinant of the (i + 1) x (i + 1) matrix obtained from T-^i by replacing the last row with 
(1, A, . . . , A*). The following simple lemma is reminiscent of the construction of classical orthogonal 
polynomials on the real line and the unit circle (see, e.g. Q) and will be useful for us in what 
follows. 



Lemma 3.4 Let X be invertible matrix with a moment sequence (3.1). For m G Z, define a 
Laurent polynomial 

R{X) = X'^vf\x) . (3.18) 

Then 



efR{X)X°'ei = {a = I + I - m - i, . . . ,1 - m) 



e'{R{X)X^-"'+^ei = (-l)MetT.^^+ 



efi?(X)X'~"-Vi = det7;^|\ . 



(3.19) 
(3.20) 
(3.21) 



Proof. -R(A) can be written as X]/3=o -^/3-^"^^^' where is equal to (— 1)*+^ times the minor 
obtained by deleting the last row and (/3 + l)st column in ( 3.17| ). Then 



ejR{X)X'^e,=J2Rp 



(3.22) 



If a is in the range specified in (3.19), the right-hand side of ( 3.22| ) becomes a determinant in which 
two of the rows coincide. Equalities ( 3.20| ), ( ^.21 ) are obtained in a similar way. □ 
Lemmas above imply the following 

Proposition 3.5 Assume that Af'^ ^ for i = 1,. . . ,n. Then polynomials defined by formulas 



A 



1 



(0) 



(3.23) 



satisfy ^ ^) 
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Proof. IfK(A) = EI =QPii^^ is a solution of [3l^, |3^ ), then by Corollary ^ and (|3l] ), we have 

i i 

eJxyi{X)ei =0 (1 - ^ e/3 < a < z - 5^ ) . (3.24) 

/3=1 ^=1 

If Af ^ / 0, Lemma |^ imphes that ( |3^ ) is a unique monic polynomial of degree i that satisfies 
( 3.24 ). Therefore it has to coincide with the unique solution of (|3.2D . This completes the proof. □ 
Now we will be able to express parameters Cj, di in terms of Toeplitz determinants ( p. 16 ). First, 
we need the following elementary lemma. 

Lemma 3.6 Let X be an arbitrary upper Hessenberg matrix and let Dq = 1 and Di (i = 1, . . . ,n) 
denote the left upper i x i principal minor of X. Then 

Di = (-l)Vi(O) , (3.25) 

where Pi{X) are polynomials defined by ( \3.^) . 

Proof. The Laplace expansion w.r.t ith column leads to the following recursion for Di: 

i 

Di = Y^i-irxaiD^.i 

or 

i 

Di + ^(-l)"-ixa^Z)a-i = , 

a=l 

which coincides with the recursion for pi{0) obtained when one sets A = in (|3.2| ). □ 
Theorem 3.7 

_A«aS _ Al^r^Aji^') 

i j— 1 i i 



Proof. First note that, by ( 3.23 ), the right-hand side of ( 3.25| ) is equal to — whereas, due to 

(2.5), the left-hand side is di • • • dj. This immediately implies the first of the formulae ( 3.26 ). 

Secondly, it is clear from ( |3.8|) , that parameters are uniquely determined by di and 7^. On 
the other hand, ( |3.8D , together with ( |3.6D , gives 



7i+i = e'[X''^+^pi{X) 



ei 



(3.27) 



Recall that, if ej+i = 1 then i/j+i = — e/3 and i — X]^=i = i + ^ ~ Z]/3t=\ ^13- Otherwise, 

i/j+i = i + 1 — Y2^f^=i = ^ + 1 — ^''fi^i e/3- Comparing ( |3.27D with equalities ( 3.2[l| ), (|3.21| ) and then 
using ( p. 16 ), one obtains 



7i+i = (-l)->+icidiii^^^-^— ^ = (_i)Ki-»+i)^ 



(di 



■di+i)^'+i 



A 



(0) 



(3.28) 
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Now, to finish the proof, it is enough to check that when one substitutes expressions ( p.26 ) into 
the right-hand side of (l3.8|) , the result agrees with ( 3.28 ). It is not hard to see that if dp,cp are 
defined by ( |3.26| ) then 



a(^/3+i) .(0) 



(3.29) 



^(^2) 



Ah) 



while diCi = j^-. Therefore, the numerator in the right-hand side of (3.8) is equal to —for ^"^^ 



(|3^ ) and (|]2|) agree. 



denominator is equal to ( ^^gy ) . For both possible values of ej, this implies that formulae ( p.8| 



□ 



One of the consequences of Lemma 3.2 is an existence of the unique sequence of Laurent poly- 
nomials ro(A), . . . , r„_i(A) of the form 



n{\) = X'^^' +Y,rip\''^' ■ (3-30) 

/3=1 

such that 

efri_i(X) =7ief (i = l,... ,n) (3.31) 
Next proposition shows that can be conveniently described by formulas similar to ( p. 23) ). 
Proposition 3.8 



A 



(0) 



'(A) 



(3.32) 



Proof. First note, that by ( ^l3|) , (HJ^) can be re-written as 



ri(A) = Ai-^?=i^'' J^a^^AV (3.33) 

/3=1 

Next, by (p.ll|) , fj+i coincides with the lowest (resp. highest) degree in rj(A) if e^+i = 1 (resp. 
e^+i = 0). Finally, due to Lemma 3.1, the claim that e^rj(X) is proportional to e^j^i is equivalent 
to a property 

elri{X)X^ei = (Z = 0, . . . , i - 1) (3.34) 
This property is satisfied by Lemma □ 
Corollary 3.9 Define a bilinear form (,) on C[A,A^"'^] hy 

(V,A^') = s,+j (i,i GZ) . 
Then functions po{\) ,.. . ,Pn-i{^)', r'o(A),... ,rn-i{X) form a hi-orthogonal system w. r. t. (,), i. 



e. 



(3.35) 



Proof. Follows immediately from Lemmas 3.1, |3.2| and Propositions 3.4, 3. 



□ 
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4 Symmetric Case 



It is natural to ask how should results of the previous section be modified, if initially one identifies 
bl with a space S of symmetric matrices rather than with 7i. In this case, a co-adjoint orbit of 
B_ through Xq € 6 is described as 



OJsym 



: n G B_} , 



(4.1) 



where A 



sym 



A>0 + (^>o)^- 



A set Mj is still defined by ( |2.3| , but the definition of Xj should be changed as follows: 

fc-i 



(4.2) 



An open dense subset Mj C Mj, that we are going to study, consists of elements of the form 

X = (1 - U^r\l + Ul)D{l + C/i)(l - U2)-\ (4.3) 

where, as before, D = diag((ii, . . . ,^2) and J7i, U2 are defined by ( 2.10| ). Matrix entries of X then 
are found to be 



l-l 



Ci ■ ■ ■ Cm-l{dl + ^ da{c 



Ij-l 



< I < m < i-i 



1 



dij_iCij_-^ ■ • • CfYi—i, 



ij-l = I < m < ij 
otherwise 



(4.4) 



A similar expression can be obtained for matrix entries of X ^. Denote by J a set of indices 
({1, . . . , n} \ /) U {n} = {/o = 1< /i < • • • < In-k = n}. Then 



ix-')i^ = (x-i) 



ml 



Cri 



L-i—i 

,1 ^ 1 

tLm. Urv- 



0, 



Ij-i < I < m < Ij 
Ij-i = I < m = Ij 



otherwise 



(4.5) 



Note that a conjugation of X by the matrix diag(lj , — l^-i) does not change values of parameters 
dj,j = 1, . . . ,n and Cj,j / i, but changes q to — q. This conjugation also does not affect values 
of the moments of X. Thus, in order to make a solution of the inverse problem below unique, we 
shall assume that all Cj are positive. 

As in the previous section, we are interested in expressing parameters Ci,di via the moment 
sequence (sj = Si{X) = ef X^ei)i^z of X. 

We start by noting that Lemma 3.2 remains literally true for matrices X of the form 4.3, which 
implies an existence of the unique sequence of Laurent polynomials ro(A), . . . ,r„_i(A) of the form 
( 3.30| ) satisfying ( 3.31 ). Formulas for functions rj(A) are similar to ( 3.32| ). 

First, we define a new collection of Toeplitz determinants 



(m) 



det r, 



(i+l+m-2 5:jj^-^e^) 



(m G Z) . 



(4.6) 
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Proposition 4.1 



Proof. We can argue exactly as in the proof of Proposition |3.8| . Note, however, that since X is 
symmetric, rj_i(X)ei = 'yiCi. This means that condition ( 3.34] ) that guarantees ( |3.31| ) has to be 
replaced by 

eJri^i{X)X''"ei = (a = 1, . . . , i - 1) 

or, equivalently, 

i-l i-l 

efri_i(X)X'ei = 0(Z = l-5^e^,... ,i-l- J^e^j) . 

/3=1 /3=1 

Function rj_i defined by ( |4.7| ) satisfies this condition for exactly the same reason that function 
djU) satisfies ( FP) . □ 



Corollary 4.2 



2^(0) 



7f = (-1)^-^^ (4i 



Proof. It follows from ( |3.7| ), ( p.31 ) and symmetricity of X that 

7,2 = efri„i(X)X^>ei . (4.9) 

Then (|]|) follows from Lemma 0, ( pl| ) and (|]]). □ 

Remark. Note that, as one would expect, in the case of the classical moment problem (I = 
{2, 3, . . . , n}, X is tri-diagonal) , a condition that ensures that the right-hand side of (4.8) is positive 
is a positive definiteness of the Hankel matrix {si+j-2)i'j=i- Similarly, in the case / = {n} which 
relevant in the study of the peakons lattice and was comprehensively studied in |2j, one comes to 
a conclusion that the Hankel matrix {s2~i-j)fj=i is positive definite. 

Let us consider now an m x m sub-matrix of X E Mj obtained by deleting (n — m) last 
rows and columns. It is clear from (|4.4|) , that Xm does not depend on parameters Cm, • • • , Cn-i, 

dm+l ) • • • ) ■ 

Lemma 4.3 Let Sa{Xm) = eiX°'ef . then 

m m 

So,{X„,) = Sa{X) (2-2^e^<a<2m + l-2^e^) . (4.10) 

/3=1 /3=1 



Proof. For / > 0, an expression for si{X) in terms of matrix entries of X reads si{X) = 
a^iQia^aiQ2 ■ ■ ■ ^ai-ii- (iil)) many of the terms in this sum are zero. Moreover, if 
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/ < 2k + 1, were k is the cardinality of I, one can find among the non-zero terms the one 

9 

were max(ai,... ,a/_i) reaches its maximum. This term is equal to JJ^?a-i,«/3' if ^ = and 

/3=1 

Q 

■^iqiq W_ -^"ip 1 i^) if ^ = 2g + 1. This implies that if m > ij and < I < 2j + 1 then the expression 

/3=1 

for si{X) involves only entries of and, therefore, si{X) = si{Xm)- Since the largest j such that 
m > ij is given by m — X]/?=i ^/3' ( 4.1C| ) is satisfied for a = 0, . . . , 2m + 1 — 2 Yl^=i ^P- 

Similarly, one can use ( fl.SD to conclude that if q is the largest index such that Iq < m then (i) 
for a,/5 < Ig, {X~^)af3 = {X~^)af3; (ii) for —2q < I < 0, the expression for si{X) contains only 
entries {X~^)ai3 with a,/3 < Ig, which means that in this case si{X) coincides with si{Xm)- To 
finish the proof, it remains to notice that q is equal to Yl^=i ~ ^- ^ 



Proposition 4.4 



det(A - X„) = (A) • (4.11) 



Proof. Let = A™ + X^^^q^ QmiA*. Then the Hamilton-Cayley theorem implies 



m—l 



{Xm) + 

(X^) = (a G Z) . 



(4.12) 



i=0 



By Lemma [4.3| , ( |4.12D remains valid if we replace Sa+i{Xm) with Sa+i = Sa+i{X) for i = 0, 
as long as 2 — 2 e^<a<m + l — 2 . This means that, after the right multi- 



plication of the matrix used in the definition ( p.l7| ) of V,^ 



{m+l-2E^=i€/3) 



(1 + jyi^^i ^f3+i,m+i), one gets a matrix of the form 



m+l 



by the unipotent matrix 



{m+l-2E™.ie;3) 
m 

1 A • • • A 



m— 1 





det(A - X„ 



and (|4.1lD follows. 

It drops out immediately from ( 4.11| ) and ( [4. 3D that 



□ 



di ■ ■ ■ d„ 



(4.13) 



Now the analogue of Theorem 3.7 for matrices ( [4.3[ ) drops out immediately from ( [4.13| ), ( [4.8| 
and (p 



Theorem 4.5 



,{0) ^(1) \2 



(1) \ ^i+i 



(0) ■ 

i-l 



3 



(1) 



(0) 

i+i. 



2) 



(1) 

i-l, 



(4.14) 
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5 Conclusion 



The results of this paper can be used to make several constructions used in the study of integrable 
lattices of Toda type more explicit. For example, in |Q we proved that for any /, J there exists 
a bi-rational Poisson map from M| to Mj, that intertwines the Toda flows on and M'j. A 
construction we gave was by induction and expressions for matrix entries of elements of M'j in 
terms of matrix entries of elements of M[ were very involved. Now we can give explicit formulas 



for this map using just (3.1) and Theorem |3.7| . Similarly, Theorem 4.5 can be used to simplify ( 
in the case of elementary orbits ) the construction of the Poisson map from the Kostant-Toda to 
symmetric Toda flows proposed in Q. 

A solution to the inverse problem can also be useful in avoiding blow-ups of the Toda flows 
by switching from one elementary orbit to another, which, in turn, may become important in 
the context of LU type algorithms for computing eigenvalues of Hessenberg matrices. A natural 
question that arises in this connection is to give an intrinsic description of all sequences (sj) from 
which an element of M'j can be restored for at least one /. 

Because of possible implications in the coding theory, it would be also interesting to study 
solvability of the inverse problem over a finite field (cf. where the connection with the coding 
theory was observed in the tri-diagonal case). 

Another direction of possible investigation that we plan to pursue in the future is an extension 
of our results from M'j to M/ (to this end, some of the genericity assumptions will have to be 
relaxed) and then to more general co-adjoint orbits. In the latter case, one has to restore and 
element of the orbit from a collection of rational functions. In the case of generic co-adjoint orbits 
this has been done in ||6|, [l^ . 
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